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Abstract 

Multiple permanent-wave trains in nonlinear systems are constructed by the asymptotic tail- 
matching method. Under some general assumptions, simple criteria for the construction are pre- 
sented. Applications to fourth-order systems and coupled nonlinear Schrodinger equations are 
discussed. 



1 Introduction 



Nonlinear wave systems have been studied for a few decades. Much progress has been made on 
integrable equations where the inverse scattering transform method can be applied [1]. For non- 
integrable equations, the general analytical treatment has been elusive so far and will likely remain 
so in the near future. A less ambitious goal, then, would be to generally study the permanent waves 
in non-integrable systems. Such waves often contain valuable information on the system's general 
solution behaviors. An interesting fact is that, in many non-integrable systems, simple permanent 
waves can be matched together and form multiple permanent-wave trains ([2], [3], [4], [5], [6], [7], 
etc.). If solitary waves exist, multiple solitary-wave trains can be constructed by a perturbation 
method proposed by Karpman and Solov'ev [8] and Gorshkov and Ostrovsky [9] (see [3]). If perma- 
nent waves with exponentially decaying and oscillating tails are present, the existence of countably 
infinite multiple permanent-wave trains has been proved for certain types of nonlinear systems by 
variational methods ([4], [5]). In this paper, if a nonlinear wave system allows permanent waves 
which exponentially approach a constant at infinity, we will construct widely-separated multiple 
permanent- wave trains by a new and general method, namely, the asymptotic tail-matching method. 
This method is stimulated in part by another matching method for non-local solitary waves (see [6] 
and [7]). Under some general assumptions, we will show that an arbitrary number of permanent 
waves can be matched together and form multiple permanent- wave trains if and only if the expo- 



nential tails of these permanent waves satisfy certain simple algebraic conditions. These conditions 
will also determine the spacings between adjacent permanent waves if such matching takes place. 
This asymptotic tail-matching method differs from Karpman et al's perturbation method in two 
major aspects. First, it can be applied directly to the matching of kink and anti-kink type perma- 
nent waves. Second, its results are explicit, simple and insightful. As applications of these general 
results, we will discuss fourth-order systems and the coupled nonlinear Schrodinger equations. For 
fourth-order systems which allow permanent waves exponentially and oscillatorily approaching a 
constant at infinity, we will show that countably infinite multiple permanent- wave trains exist and 
can be readily constructed. Thus the results in [4] and [5] are reproduced. For the coupled nonlin- 
ear Schrodinger equations, we will show that countably infinite multiple solitary- wave trains can 
be constructed in a large portion of the parameter space. Numerical results will also be presented 
and compared with the theoretical predictions when appropriate. 



2 Construction of multiple permanent-wave trains 

We consider a general nonlinear wave system 



where U is the unknown vector variable, and F is a nonlinear vector function. Suppose it al- 
lows permanent waves of certain form which, when substituted into Eq. (2.1), reduces it into an 
autonomous complex system of first-order nonlinear ordinary differential equations 



where $(x) is a n-component vector variable. If Eq. (2.2) has permanent wave solutions which 
exponentially approach a constant at infinity, then we next will develop a new method to determine 
if those permanent waves can be matched together and form widely-separated multiple permanent- 
wave trains or not. The idea is to perturb each permanent wave such that the exponential tails of 
each perturbed wave match those of the adjacent permanent waves. We first discuss the matching 
of solitary waves, followed by that of general permanent waves. 

2.1 Solitary- wave trains 

Suppose Eq. (2.2) allows solitary waves <&(x) which exponentially decay to zero as |x| — > oo. We 
make the following general assumptions: 

Al. the eigenvalues of the constant (Jacobian) matrix VG(0) all have non-zero real parts; 

A2. for any solitary wave &(x), the linear behavior dominates at infinity, i.e., as x — > oo or — oo, 
$(x) approaches a solution of the linear equation 



F(U, D X: D t ) = 0, 



(2.1) 



d<S>/dx = G($) 



(2.2) 



d$/dx = VG(0)$; 



(2.3) 
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A3, for any solitary wave &(x), the linearized equation of (2.2) around <&(x) 

d$/dx = VG($)$ (2.4) 

and its adjoint equation 

-d^/dx = VG* T ($)^ (2.5) 

each have a single linearly independent localized solution. Here "T" represents the transpose 
and "*" the complex conjugate. 



Remark: Since Eq. (2.2) is autonomous, any spatial translation of <&(x) is still (2.2)'s solution. 
Therefore Eq. (2.4) always has a nontrivial localized solution d<&{x)/dx. The requirement for 
Eq. (2.4) is just that d<&{x)/dx is its single linearly independent localized solution. This can be 
guaranteed if &(x) is isolated in ^(R, R?) up to spatial translations (the so-called non-degenerency 
condition in some of the literature. See [4]). 

We also introduce the following notations. In view of assumption Al, let us denote VG(0)'s 
eigenvalues as Ai, A2, . . . , A n where 

Re(Ai) < Re(A 2 ) < . . . < Re(A s ) < < Re(A s+ i) < Re(A s+2 ) < . . . < Re(A n ). (2.6) 

Suppose for an eigenvalue A, VG(0) has a chain of eigenvector and generalized eigenvectors V{ (i = 
such that 

(VG(0) - Ai>i = 0, (2.7) 
(VG(0) - M)v i+1 = Vi , i = l,...,l-l. (2.8) 
Define the polynomial functions £i(x) (i = 1, . . . , I) as 

-1 

(i-Ty: 

then 

£+10*0 = 60*0, (2-io) 

and {£i(x)e Xx , i = 1, . . . , 1} form a chain of linearly independent solutions of Eq. (2.3). According 
to the theory of linear differential equations with constant coefficients, we can find such chains 
of solutions which together form a fundamental set of solutions of Eq. (2.3). Thus according to 
assumption A2, we have 

HX) \ E?= s +i Cittx)***, x -00, (2 - U) 

where c, (i = 1, . . . , n) are complex constants. We point out that in the special case where VG(0) 
has n linearly independent eigenvectors, i = 1, . . . , n} are just those constant eigenvectors. The 
fundamental matrix of the adjoint equation (2.5) at infinity is 

[ m e-^ x m e-^ x ... r! n e- x >], (2.12) 



£i(x) = Vi + xvi-i + . . . + X Tr: vi, i = l,..., I, (2.9) 
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where 

lm m ••• vn] = {[ti 6 ■■■ CrT 1 }^ (2.13) 

Note that for 1 < i,j < n, 

e*(x)-nj(*) = { J; -/j; (2.i4) 

Thus for the single linearly independent localized solution ^(x) of Eq. (2.5), we have 
where di (i = 1, . . . , n) are complex constants. 

Now suppose {$(!), $( 2 ) $W } are iV solitary waves of Eq. (2.2) with 

3W (x) _^ { £U ^to^' * ^ °°> (2.16) 

For each <£>( fc ), the single linearly independent localized solution of the adjoint equation 

-dV<V/dx = VG* T (<5>W)¥^ (2.17) 
has the following asymptotic behavior at infinity: 

• w <*> - { ^- 'Of/t A 1 " (2 ' 18> 

Consider a new solitary wave which looks like a superposition of the above N solitary waves {$( fc )} 
widely separated, with the fc-th wave $( fc ) located at x = x& (k = 1, 2, ... , TV). Let 

xi < x 2 < . . . < xat, (2.19) 

and denote 

A fc = x fc+1 - x k (» 1), fc = l,2,...,JV-l. (2.20) 

We will call this new solitary wave as a TV-pulse wavetrain. It can be constructed explicitly by the 
following theorem. 

Theorem 1 Under the assumptions Al, A2, A3 and the above notations, the N solitary waves 
{(E^ 1 ), . . . , cfrW} can match each other and form a widely-separated N-pulse wavetrain if and only 
if the spacings (S> 1) (k = 1, . . . , N — 1) asymptotically satisfy the following N conditions 

J2 cfdf*e- x ^ = 0, (2.21a) 

j=8+l 

^cf-^df^e^^-i = jr cf +1) df ] *e-^ A \ (2<k<N-l), (2.21b) 

j=l j=s+l 

^ c (^-i) d W* e W-^=0. (2.21c) 
The relative errors in Eqs. (2.21) are exponentially small with the spacings. 
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We will prove this theorem by the asymptotic tail-matching method to be developed next. 

Proof: Suppose such a iV-pulse wavetrain <3>(x) exists. Then around the fc-th wave (2 < k < N — 1), 
the solution is 

= - x k ) + <l (/c) (x - x fc ), (2.22) 

where <C 1. The linearized equation for is 

d& k) {x)/dx = VG($ (fc) (s)) l> (fc) (x). (2.23) 

According to assumption A3, Eq. (2.23) has a single linearly independent localized solution which 
is d& k \x)/dx. From Eq. (2.16) we get 



d ^\x)/dx fl 4 k Hzfc) + x M x )) x ^™- ( 2 - 24 ) 
=i 



Clearly not all the c^'s (i = 1, . . . , s) are equal to zero. Without loss of generality, we assume that 
c[ k ^ / 0. Then we denote the other n — 1 solutions of Eq. (2.23) as (j = 2, . . . , n). We require 
that 

$f \x) — ► (j(x)e XjX , x^oo. (2.25) 

As x — ► — oo, we generally have 

n 

$f) {x) _^J24ki(x)e XiX , i = 2,...,n, (2.26) 
i=i 

where djf are constants. Since c[ k ^ / 0, these n solutions {d&^/dx, ■ ■ ■ , are linearly 

independent at x equal to infinity, and they form a fundamental set of solutions of Eq. (2.23). 
Thus the general solution for &( k \x — x k ) is 



& k \x - x k ) = ^——(x -x k ) + J2 hM\x - x k ), (2.27) 

dx j=2 



where hj (j = 1, . . . , n) are constants. The first term in (2.27) can be absorbed into <&( k \x — x k ) 
and cause a position shift to it. By normalization we make h\ = 0. When x k <C x <C x k+ ±, dropping 
the exponentially small terms, we get 

n 

$(*>(z-z fc )— ► Yl hjij{x - x k )e x ^ x - Xk \ (2.28) 

j=8 + l 

Similarly, when x k -\ <C x <C x&, we have 

*(*)( x _ Xfc ) _^ ^EfllfhjJ^s - x k )e x ^- x "\ (2.29) 

The key idea in the asymptotic tail-matching method is that, in order for the matching to occur, we 
need to require that in the region x k <C x <C x k +i, ^ k \x — x k ys exponentially growing terms match 
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the left tail of the right-hand wave cE>( fc+1 )(x — ajfc+i); in the region x^-i <C x <C x&, <l>( fc )(x — Xfc)'s 
exponentially decaying terms match the right tail of the left-hand wave <I>( fe_1 )(x — Xk-i)- In the 
region x& <C x <C Xfc+i, this requirement is 



E ^(x-x fe )e^^)= £ cf +1) ^(x-x fc+1 )e^^+^; (2.30) 

j=s+l j=s+l 



and in the region x^-i <C x <C x^, it is 

s n s 

E(E - x fc )e^(— *) = ]T cf - 1} ^(x - x^Oe^'^-i). (2.31) 

j=l i=2 i=l 

We now need to select the constants hj (j = 2, . . . , n) and spacings (A; = 1, . . . , N — 1) so that 
the above two conditions are satisfied. 

First consider condition (2.30). Recall that functions {£j(x)} are of the form (2.9). If for A = 
X m (s + 1 < m < n), the chain of such functions is {£ m , . . . , £ m +/_i}, where £ m is a constant vector 
and 

^H-i(z) =&(x), j = m,...,m + l-2. (2.32) 
Then we select {h m , . . . , h m+ j-i} from the equation 

m+l—l m+l—1 

E h^(x-x k )e x ^-^= £ cf +1) ^(x-x fc+1 )e A ^+i). (2.33) 

jr'=m j=m 

The right hand side of this equation is 

m+l-l 



E 4 fe+1) ^(x-x fc -A fc )e^)e-^ 



(fc+i), 

= E # +1) (E^^-**)>«^-^- Ai » 

m+/-l j / a 

= E E4 fc+1) T^r^(--^) eA( ™ )e - AAfc 

m+l—1 m+l—1 



= E ( E V^yH i ) e "^( a -^ Mi) - ( 2 - 34 ) 

i=m j=i \J ' 

Now we choose hi (i = m, . . . , m + I — 1) to be 

^ = (E f +1) )e- AAfc , i = m,...,m + Z-l, (2.35) 



3=1 



then Eq. (2.33) is valid. Repeating this procedure for the other chains of {£j(x)} functions in the 
form (2.9), we can successfully select hi (i = s + 1, . . . , n) so that condition (2.30) is satisfied. 
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Next consider condition (2.31). Similar analysis shows that we can reduce its right hand side to 

E c f~ 1] ti( x - Xk-i)e X ^ x - x ^ = E aje^-^jix - x k )e x ^~^\ (2.36) 

3=1 3=1 

where ctj (j = l,...,s) are constants and determined by c^ -1 ^ (j = l,...,s) and Ak-i- Then 
condition (2.31) becomes 

j^ a f) hi = aje ^ k -i- £ a f hh j = l,...,s. (2.37) 

i=2 i=s+l 

This is a linear system of s equations for s — 1 unknowns hi (i = 2, . . . , s). We now show that 
the matrix (o^) sx (s-i) on the left side of Eq. (2.37) has rank s — 1. Consider the solution of Eq. 
(2.23) 

T(s) = Eft-*f (*). (2-38) 

J'=2 

where pj (j = 2, . . . , s) are constants. Dropping exponentially small terms we get 

I x — ^ oo 

T(s)—> {' r*)' t ,^, x „ _ (2-39) 



E?=i(Ef=2< ; Pi)^Wc AiX . 



x — > — oo. 



According to assumption A3, the only localized solution of Eq. (2.23) is d& k \x)/dx. Moreover, 
c[ k ^ in (2.24) is non-zero. Thus T(x) can not be a localized solution. In other words, the linear 
system of equations 

s 

E a lJ ) Pi = °. J = l,...,s, (2.40) 

i=2 

has no non-trivial solutions for pi (i = l,...,s). Therefore the matrix (a^)sx(s-i) has rank 
s — 1. Without loss of generality, we assume that the last (s — 1) rows of the matrix are linearly 
independent. Then the linear system 

s n 

£ afhi = aje^-i - £ afh t , j = 2,...,s, (2.41) 

i=2 i=s+l 

has a unique solution for hi (i = 2, . . . , s). With hi (i = 2, . . . , n) given by (2.35) and (2.41), the 
only matching condition left to be satisfied now is 

£ 4?hi = a ie ^- - £ o.?^, (2-42) 

j=2 i=s+l 

which will determine the spacings of this iV-pulse wavetrain. Since the matrix (a^) sx (n-i) is n °t 
readily available, to determine the spacings from Eq. (2.42) is difficult. But this can be easily done 
with the aid of the solution of the adjoint equation (2.17). With hi (i = 2, . . . , n) given by 
(2.35) and (2.41), it is easy to show that Eqs. (2.28) and (2.29) become 



¥ k \x-x k )^ cf +1 kj(x - x k+1 )e x ^ x - x ^ , x fe «x«x fc+1 , (2.43) 

j=8+l 
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and 



J'=2 

(2.44) 



where w is a constant. Condition (2.42) is equivalent to 

w = cf~ 1] . (2.45) 

For Xk-\ <C y\ <C x& and < 1/2 < iCfe+i, we have 

= j^ 2 {dl>( fc )(x-x fe )/dx-VG($( fc )(x-x fe ))$W(x-x fc )}-^W*(x-x fc )dx 
= $«(x-x fe ))-^ fc )*(x-x fc )|^ 

+ Pp^-^.j-^Wfj-Xij/dx- VG* T ($( fc )(x-x fc ))^( fc )(x-x fc )}*(ix 1 " j 
= m\x-x k ))-^ k >{x-x k )\y\. 

For Afe_i » 1 and A*, 3> 1, asymptotically we get 

w d^V 1 **" 1 + 2 C f-V d f*e x ^ = J2 cf +1) df*e- x ^K (2.47) 

j=2 j=8+l 

Condition (2.45) is satisfied if and only if Eq. (2.21b) is valid. For the first and last waves in this 
TV-pulse wavetrain, the analysis is simpler, and we get Eqs. (2.21a,c) for matching. In summary, 
the N pulses -f^ 1 -*, . . . , can be matched and form a iV-pulse wavetrain if and only if the 

spacings Afc(3> 1) (k = 1, . . . , N — 1) asymptotically satisfy the TV conditions (2.21). 

Now we discuss the accuracy of the above results. Error is created mainly by the matching require- 
ments (see Eqs. (2.30) and (2.31)) and the negligence of nonlinear terms in Eq. (2.23). First we 
discuss the error in the matching requirements. Let us reconsider the solution (2.22) around the 
k-th. wave. When x k <C x <C x k+ i, beside the exponentially decaying terms in &( k \x — x k ), there 
are also such terms in & k \x — x k ) (see Eq. (2.27)). The combined exponentially decaying tails are 

cf\i{x - x k )e x ^~^ + £(cf + hj)ij{x - x k )e x &- x >>\ (2.48) 

i=2 

(k) (k) 

Thus Eqs. (2.21) would be more accurate if the Cj values are replaced by Cj + hj. Recall that 
hj (j = 2, . . . ,s) are determined from Eq. (2.41), so they are exponentially small for large A&_i 
and A k . As a result, the negligence of tail contribution from &( k \x — x k ) only causes exponentially 
small relative errors in Eqs. (2.21). Simple reasoning also shows that the exclusion of nonlinear 
terms in Eq. (2.4) also causes only exponentially small relative errors in (2.21). The proof of 
theorem 1 is now completed. It should be pointed out that, if the eigenvalues {Aj, i = 1,... ,s} 
or {Aj, i = s + 1, . . . , n} are real-valued and close to each other, those exponentially small relative 
errors in Eqs. (2.21) may become significant. In such cases, caution is needed in interpreting the 
results from (2.21) 
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2.2 General permanent- wave trains 



The results in the previous section can be readily extended to the matching of permanent waves 
which exponentially approach a complex constant at infinity. Suppose such permanent waves exist 
in Eq. (2.2), then we make the following general assumptions: for any permanent wave <J?(x) where 

*(*) — { J 2 ' X ^ °°' (2.49) 
I 01, x — > -co, 

Bl. the eigenvalues of the constant (Jacobian) matrices VG(oi) and VG(&2) ah have non-zero 
real parts, and the number of VG(oi)'s eigenvalues with negative real parts is equal to that 
of VG(02)'s eigenvalues with negative real parts; 

B2. the linear behavior dominates at infinity, i.e., as x — ► — oo and oo, $(x) approaches a solution 
of the linear equations 

dkjdx = VG(oi)$, (2.50) 

and 

d^/dx = VG(6 2 )4>, (2.51) 

respectively; 
B3. the linearized equation of (2.2) around <3?(x) 

d$/dx = VG($)$ (2.52) 

and its adjoint equation 

-d^/dx = VG* T ($)V (2.53) 
each have a single linearly independent localized solution. 

Now suppose {^ 1 \^ 2 \ . . . , $( iV )} are N permanent waves with 

r ilk) 

<S> (k) (x)-^{ Ik) (2-54) 

I, »! , X — > — OO, 

where 1 < k < N. If they are to be matched and form a widely-separated A^-permanent-wave train, 
we need to require that 

b[ k) = b ( 2 k ~ 1] , 2<k<N. (2.55) 

One consequence is that all the matrices VG(b[ k ^) and VG(b^) (1 < k < N) have the same 
number of eigenvalues with negative real parts, which we denote as s. We introduce the following 
notations. Denote VG(i)f')'s n eigenvalues as X 2 k \ • • • , A^ with 

Re(Af } ) < Re(A^) < . . . < Re(A^) < < Re(A^ ) 1 ) < Re(A^ 2 ) < . . . < Re(A^), (2.56) 

and VG(o 2 fc) )'s as A{ k) ,A ( 2 k \. . . , K {k) with 

Re(AS fc) ) < Re(A^) < . . . < Re(A<*>) < < Re(A^) < Re(A^ 2 ) < . . . < Re(A^). (2.57) 
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The fundamental sets of solutions of the linear equations 

dkjdx = VG(b[ k) )$ (2.58) 

and 

d^/dx = VG(4 fe) )4> (2.59) 

are respectively {^ k \x)e x i )x , i = l,...,n} and {6\ k \x)e A i )x , i = 1, ...,n} which consist of 
chains of linearly independent solutions of Eqs. (2.58) and (2.59) as defined before (see Eq. (2.9)). 
The fundamental matrices of the adjoint equations of (2.58) and (2.59) are then 



and 

with 

and 

Note that 



[^e-W* vi k) e- X > x ... vi k) e- X " x ] (2.60) 

[C W e -Ai fc) *x C (*) e -4«-x ... C^e-^'l, (2.61) 

# ••• ^{K? 5 ••• ^l- 1 }" (2-62) 

[ C { fc) cf ... & = {[e{ k) 0? ... eW]- 1 }". (2.63) 



(*) = <*) = (*) = . = x> _ _ _ >n (2 _ 64) 



in view of (2.55). In those notations, we have 



$(*)( X ) _> 6 2 + Ei=i q J (x)e . x - oo, ( } 



according to assumption B2. For the single linearly independent localized solution vl/W of the 
adjoint equation 

-dftW/dx = VG* T (<S> {k) )¥ k \ (2.66) 

we have 

— { El=1 d ^i)$ e ~ X ^k X " (2-67) 

Here and df^ (1 < i < n) are complex constants. 

Now consider a widely-separated permanent- wave train matched by the above N permanent waves 
$(2) ) . . . > Assume that the fc-th wave $ (fc) is located at x = x k (fc = 1,2,... , AT), and 

Afc is as defined in Eq. (2.20), then we have the following result. 

Theorem 2 Under the assumptions Bl, B2, B3, (2.55) and the above notations, the N permanent 
waves can match each other and form a widely- separated N -permanent-wave train 
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if and only if the spacings (^> 1) (k = 1, . . . , N — 1) asymptotically satisfy the following N 
conditions 

£ cf >#>V A ? )Al = 0, (2.68a) 

j=s+l 

J2 cf~ l) df*e X t )A ^ = jr cf +1) df*e- K T A \ (2<k<N-l), (2.68b) 

j=l j=s+l 

^cf-'U^e^^-^O. (2.68c) 

3=1 

The relative errors in Eqs. (2.68) are exponentially small with the spacings. 

The proof for this theorem is similar to that for theorem 1, and is thus omitted here. 

Remark: In applying theorem 2 to a given nonlinear wave system, the major difficulty is the 
determination of the coefficients {d^} in the localized solution ^( k \x) of the adjoint equation 
(2.66). In general, this has to be done numerically. But in many cases, Eq. (2.52) can be cast into 
a self-adjoint system (see [2], [4] and [5]). Then ^( k \x) and its coefficients {d^} can be readily 
obtained from d&^/dx, and the verification of conditions (2.68) can proceed. 

3 Applications 

3.1 Fourth-order systems 

The permanent waves in many nonlinear wave problems are governed by fourth-order systems (2.2) 
(see [2], [4] and [5]). In this section, we apply theorems 1 and 2 to certain classes of such systems. 
In particular, we will establish the existence of countably infinite multiple permanent- wave trains 
under some general assumptions. 

We first consider the matching of identical permanent waves in a fourth-order system (2.2). Suppose 
<3?(x) is a permanent wave in (2.2) where $(x) — > b as \x\ — > oo, and the assumptions Bl, B2 and 
B3 are satisfied. Moreover, we suppose the eigenvalues of VG(b) are ±Ai and ±A2, where Ai / A2 
and Re(Aj) > (i = 1,2). Then corresponding to the four distinct eigenvalues — Ai,— A2,Ai and 
A2, VG(b) has four linearly independent eigenvectors £j (i = 1, . . . ,4). If we denote 

[mmmv4] = mi^^U]- 1 y T , (3.1) 

then 

$(x) - 

and 

*(x) 



b + ciCie-^ + c 2 &e- X2X , 
b + c 3 £ 3 e Aia; + c±Ue Xrx , 

j dir]ie x i x + d2r/2e x 2 x , 
I d 3 r] 3 e- x i x + d 4 T]4e- X 2 x , 



x 

X 



00, 



-00, 



X 

x ■ 



— 00 
00. 



(3.2) 
(3.3) 



11 



For some fourth-order problems, Eq. (2.52) can be cast into a self-adjoint system, and one has 
either Ai and A2 real-valued with (di d%) oc (03 c±) and (d 2 d^) oc (04 c 2 ), or Ai and A2 complex- 
valued with A2 = A*, [di ds) oc (04 C2) and ^4) oc (03 ci). In such cases, conditions (2.68) for 
the matching of N identical permanent waves {&(x), . . . , <&(#)} simply become 

Cl ^e- AlAfc + c 2 d* 2 e- X2Ak = 0, fc= l,...,iV- 1. (3.4) 

In the second case, if furthermore (2.2) is a real system, then A2 = A*, c 2 = c*, cfo = d*, and Eq. 
(3.4) becomes 

Cl ^e- iIm ( Al ) Afc +c*d ie iIm ( Al ) Afc =0, k = l,...,N-l. (3.5) 
The spacings Aj, can then be easily obtained from (3.5) as 

A fc = (axg(cidJ) + |+m fc 7r)/Im(Ai), k = 1, . . . , N - 1, (3.6) 

where is any non-negative integer. Note that in this case, the exponentially small relative errors 
in (2.68) make little difference, especially when is large. Thus we conclude that an arbitrary 
number of identical permanent waves <&{x) can be matched together and form multiple permanent- 
wave trains, whose spacings are given asymptotically by Eq. (3.6). Clearly a countably infinite 
number of such wavetrains can be formed. In the paper by Buffoni and Sere [4], they proved 
the existence of countably infinite multi-pulse permanent wave solutions for a class of coupled- 
nonlinear- Schrodinger- type equations. When those equations are cast into a fourth-order system 
of the two variables and their first derivatives, it is easy to check that they fall into the above 
category. Thus their result is a special case of ours. But differences also exist between their result 
and ours. In their result, mj, in Eq. (3.6) is an even integer; while in ours, it is any integer. This 
means that we identified twice as many solitary-wave trains as they did. For fourth-order systems 
where $(cc) — b and ^f(x) element- wise are both even or odd in x, or one of them is even (odd) and 
the other one odd (even), then (03 C4) = ±(ci c 2 ), and (£3 £4) is row-wise equal to or opposite of 
(£1 £2)- It is easy to show from (3.1) that (7/3 774) is also row- wise equal to or opposite of (771 772) 
and (g?3 di) = ±(di d 2 ). Thus the N matching conditions (2.68) also reduce to (3.4). If further 
more, A2 = A*, then we will find countably infinite multiple permanent- wave trains whose spacings 
are given by (3.6). 

Next we consider the matching of different permanent waves in a fourth-order system (2.2). Suppose 
G is an odd function of <I>, i.e., 

= -G($), (3.7) 
and $(x) is a permanent wave in Eq. (2.2) with 

•(*) — f\ (3-8) 
1 —0, x — > —00, 

then — $>(x) is also a permanent wave in (2.2). It is easy to show from (3.7) that VG(— <!>) = VG( ( I > ), 
thus VG(— b) = VG(b). Beside the assumptions Bl, B2 and B3, we also assume that VG(6)'s four 
eigenvalues are ±Ai and ±A2 with Ai / A2 and Re(Aj) > (i = 1,2). Suppose the eigenvectors 
corresponding to — Ai, — A2, Ai and A2 are denoted as = 1, . . . , 4), then we have 

*( X ) _^ / b + dtie- XlX + c 2 i 2 e~^\ x ^ cx), 

{ ' \ -b + c 3 £ 3 e Al2: + c A Ue X2X , x -»• -00. 1 ' 
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For the localized solution of the adjoint equation (2.53), we have 

V(x) — ► { d 1 r ]1 e x l x + d 2 m^ x _ K x -► -oo, 

1 ^3%e x ^ x + d^e X 2 X , x — > oo, 

where r/j (i = 1, . . . , 4) are given by (3.1). For those equations (2.2) where Eq. (2.52) can be cast 
into a self-adjoint system and one has either [d\ d 3 ) oc (03 c\) and {d 2 d±) oc (04 C2) with Ai and 
A2 real, or (d\ d$) oc (04 C2) and {d 2 ^4) oc (C3 ci) with A2 = A*, conditions (2.68) for the matching 
of permanent waves {$>, — 3>, <£, — <I>, . . .} or {— <I>, 3>, — <3>, <£, . . .} will also reduce to (3.4). When 
A2 = A*, if furthermore (2.2) is a real system, then we can show as before that such matchings are 
always possible and the spacings are given by Eq. (3.6). An infinite number of such wavetrains 
will be obtained. We point out that the fourth-order systems studied by Kalies and VanderVorst 
[5] falls into this category and is thus a special case of the above results. Here again we identified 
twice as many permanent- wave trains as they did since rrik in Eq. (3.6) needs to be an even integer 
in their result. 



3.2 Coupled nonlinear Schrodinger equations 

The coupled nonlinear Schrodinger equations govern the evolution of two interacting wave packets 
in nonlinear and dispersive physical systems [10]. These equations are particularly important in 
nonlinear optics as they govern the pulse propagation in birefringent nonlinear optical fibers [11]. 
In recent years, the experimental design of high-speed optical-soliton-based telecommunication 
systems stimulated great interest in these equations, and much work has been done on them (see 
[12] and the references therein). In particular, simple and multi-pulse solitary waves in these 
equations have been found and classified in [2]. In this section, we study the multiple permanent- 
wave trains in these equations. We primarily discuss the focusing case where solitary waves exist. 
In the end of this section, we comment on the defocusing case where dark solitons arise. 



The solitary waves in coupled nonlinear Schrodinger equations (focusing case) are governed by the 
following set of equations 

rxxx ~ ri + {r\ + /?rf )n = 0, (3.11a) 



r 2xx ~ u 2 r 2 + {r\ + (3r\)r 2 = 0, 



(3.11b) 



where r\ and r 2 approach zero as x goes to infinity, and (3 and to are positive parameters. To apply 
theorem 1 to these equations, we first rewrite them as the following first order system 



where 
and 



dU/dx = G(U), 



U = (ui,u 2 , «3,^4) T = (ri,r lx ,r 2 ,r 2x ) T , 



G(U) 



( u 2 \ 

ui - {u\+f3ul)ui 
U4 

\ u 2 u 3 - (uj + puj)u 3 J 



(3.12) 
(3.13) 

(3.14) 
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It is easy to check that the above system satisfies the assumptions Al, A2 and A3 when (3^1. 
Thus in the following we assume that (3^1. The eigenvalues of the matrix VG(0) are —1, —oo, 1 
and u, and the corresponding eigenvectors are 



(3.15) 



t 1 ^ 




f ^ 








( °\ 


-1 









i 










l 





l 


V o ) 




V J 











For a solitary wave (ri,r 2 ) of Eqs. (3.11) with 



and 



we have 



n(x) 

r 2 (x) - 



c 3 e 2 



c 2 e 
c 4 e u 



X — > oo, 
x — ► — oo, 



X 



oo, 
— oo, 

x — 
x — 



oo, 



-oo. 



(3.16) 
(3.17) 
(3.18) 
(3.19) 



VG(U) = 



\ 


1 



(3.20) 



The linearized equation of (3.12) around a solitary wave U(x) is 

dU/dx = VG(U)U, 

where U = (Hi, u 2 , u 3 , ua) t , and 

/ 1 

1 - 3it? - I3ul -2/3uin 3 


\ -2f3um 3 1 - 3u§ - (3u( / 

The single localized solution of the above equations is dU/dx. If «2 and U4 in (3.19) are eliminated 
in favor of u\ and U3, then the linear system for u\ and its are self-adjoint. The adjoint equation 
of (3.19) is 

-dV/dx = VG T {U)V (3.21) 

with V = {v\,V2,vz,V4). It is easy to see that when v\ and v 3 are eliminated from (3.21), the 
equations for v 2 and V4 are the same as those for u\ and u 3 . Thus the single localized solution of 
Eqs. (3.21) is 

V = {—Ui xx , Ui x , — U 3xx , U3 X ) T . 

At infinity, 

y, x f diriie x + d 2 me UJX , % 

* [ d 3 r] 3 e~ x + d^e - ^, x - 

where {rji, i = 1, . . . , 4} are obtained from Eq. (3.1) as 



-00, 



00, 



(3.22) 
(3.23) 



/ \ \ 



m 





V J 



( \ 




V 



1 

2 

j_ 

2lu 



1 
2 


V J 



( \ 




1 

2 

V 2w / 



(3.24) 
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and 

d\ = — 2c3, CZ3 = — 2ci, CZ2 = — 2o; 2 C4, c?4 = — 2uj 2 C2- (3.25) 
Now we consider the matching of N solitary waves {(r^Kr^), k = 1, . . . , iV} where 

r<%) — H| e x ~ X ' (3.26) 
Co 'e , x — > — oo, 



'5 



and 

C4 e , x — ► — 00. 
In view of (3.25), the matching conditions (2.21) become 

c (k) c (k+i) e -A k + ^ c {k) c (k + i) e ^ Ak = Qj fe = l ; . . . ; j\T — 1. (3.28) 

Interestingly Eq. (3.28) indicates that the iV solitary waves { (rf^ , ) , k = 1,...,N} can be 
matched if and only if all the adjacent solitary waves can. Thus the matching of solitary waves 
in Eqs. (3.11) is a "local" phenomenon. This fact would make the construction of those multiple 
solitary-wave trains much easier. In what follows, we discuss the matching of some special types of 
solitary waves. 

First we consider the matching of wave and daughter wave solutions. In such solutions, either 
T2 <C r\ or n <^T2- Without loss of generality, we assume that r2 <C r± . These solutions exist near 
the curves 

u = (y/l + 8j3-l)/2-m (3.29) 



in the (uj,(3) parameter plane [2]. Here m is a non-negative integer and m < (y/l + 8/3 — l)/2. In 
these solutions, n is symmetric; r2 is symmetric for even values of m and anti-symmetric for odd 
values of m. Suppose (fi,r2) is such a solution, then 

fi — ► cie~ |ie| , |x| -► 00, (3.30a) 

and 

/ C2e _a;x , x -> 00, . . 

r2 ^((-irc 2 e- x^-co. (3 ' 30b) 

Here C2 <C 1. Notice that if r^(x) (i=l or 2) is a solution of Eqs. (3.11), so is — rj(x). Without 
loss of generality, we require that q > (i = 1,2). Now we consider ./V wave and daughter wave 
solutions { (r ^ , ) , fc = 1, . . . , N} where 

r[ k \x) = q[ k ^ri(x), r^{x) = (^^(x), k = 1, . . . , N, (3.31) 

and qf^ = ±1 (i = 1,2). The matching condition (3.28) for these solitary waves are simply 

q{ k \{ k+1) c 2 1 e- A * + {-l) m q {k) q {k+l) u 2 cle-^ =0, k=l,...,N-l, (3.32) 

i.e. 

(fe) (fe+i) 2 2 

e -(i-)A fc = ( _ ir+ i^2 ^ fc = 1) ... )7V _ L (3.33) 
q\ >q\ c i 



15 



For these conditions to be satisfied, we need to require that uj < 1 and 

(-l) m+1 q{ k) q[ k+1) qi k) qi k+1) = 1, k=l,...,N-l. (3.34) 

Suppose (q[ k \q2 k ^) is fixed, then condition (3.34) shows that (q[ k+1 \ q^ 1 ^) can take two sets of 
values. In other words, there are two possible types of matching. Thus these ./V wave and daughter 
wave solutions can form 2^ topologically distinct solitary- wave trains. Since N is arbitrary, count- 
ably infinite multiple-pulse solitary waves will be formed. The spacings between adjacent waves in 
those wavetrains are 

At= >J^Ml, (3.35) 

UJ — 1 

which are the same throughout an entire wavetrain. As u approaches the wave and daughter wave 
boundary (yj\ + 8(3 — l)/2 — m, c\ approaches 2\/2, C2 approaches 0, and thus approaches 
infinity. The above theoretical results can be checked numerically. We first select ((3, uj) to be (2/3, 
0.85) which is close to the curve (3.29) with m equal to zero. With these parameter values, it is 
easy to find numerically that c\ and C2 as in Eq. (3.30) are equal to 2.6592 and 1.1744 respectively. 
Eqs. (3.34) and (3.35) then predict that the two wave and daughter waves (r\, r-i) and (— r\, r 2 ) can 
be matched with the spacing approximately equal to 13.0635. This is indeed the case. Numerically 
we found this exact two-pulse solitary wave and plotted it in Fig. la. The exact spacing (measured 
as the distance between the two extrema in r±) is 13.064, which is very close to the theoretical 
prediction. Next we select ((3,uj) to be (2, 0.6) which is close to the curve (3.29) with m = 1. In 
this case, we numerically found that c\ and c 2 in (3.30) are equal to 3.0386 and 0.6041. Then we 
predict from (3.34) and (3.35) that (r±, r 2 ) and itself can be matched with the spacing approximately 
equal to 10.6308. Indeed, that exact two-pulse solution was numerically found and plotted in Fig. 
lb. The exact spacing is 10.40, close to the predicted value. The predictions on other types of 
matchings were also verified with good accuracy. We point out that each multiple-pulse solitary 
wave will generate a family of solitary waves as the parameter pair (uj,(3) moves away from the 
curves (3.29). Therefore countably infinite families of solitary waves will be generated near those 
curves. 

Next we discuss mixed matchings between wave and daughter wave solutions and other types of 
solitary waves. When (uj, @) is near the curve (3.29) with m = 0, beside the wave and daughter wave 
solutions, another type of solitary waves (belonging to family D2) also exist [2]. Suppose (f : i,r 2 ) 
is a wave and daughter wave solution whose large x behavior is given by (3.30) (with m = 0), and 
(f\,f 2 ) is a solitary wave with 

f\ — ► Qie - ' 11 ', \x\ — > 00, (3.36a) 

f2 — ► 02 sgn(x) e~^ x \ \x\ — > 00, (3.36b) 

and a, > (i = 1, 2). Consider the mixed matching of the solitary waves (q\f\, ^2^2) and (q%r\, 94^2) 
where % (i = 1, . . . , 4) are either 1 or — 1 . The matching condition is 

gig 3 ciaie _A - q 2 qiUJ 2 C2a2e~ u)A = 0, (3.37) 

or 

e -(i-.)A = qjm^c^ (3 _ 38) 
9193 
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where A is the spacing. This condition can be satisfied if and only if co < 1 and the sign of (j , i(j , 2 ( 2 , 3Q , 4 
is equal to 1. As an example, we choose ((3,u>) as (2/3, 0.78). Then it is easy to find that c\, 
C2, a\ and a 2 are 2.7967, 0.5210, 7.8105 and 8.4171 respectively. The above results predict that 
(fi,r2) and (— fi, —f 2 ) can match each other and form a new two-pulse solitary wave. This was 
verified numerically. The exact matched solution is plotted in Fig. 2 with the spacing 10.26, while 
the predicted value for the spacing is 9.5571. Mixed matching between many copies of (n,/^) and 
{fi,f 2 ) can be similarly analysed. Once again, countably infinite multiple-pulse solitary waves will 
be formed by these mixed matchings. 

Lastly we discuss the matching of solitary waves near uj = 1. In this case, single- hump solitary 
waves with r\ ps r 2 are present. Suppose {f\,f 2 ) is such a solution with 

fx — ► cie"l x l, |x| -► 00, (3.39a) 

f 2 — ► C2e~ w|x| , \x\ -► 00, (3.39b) 

then ci ~ C2- If we consider the matching of these solitary waves {{q^fi, q^r 2 )} where qf^ = 
±1 (i = 1,2), the matching condition would again be Eq. (3.32) (with m = 0). But here since 
VG(0)'s eigenvalues 1 and uj are close, the exponentially small relative errors in (2.21) and (3.32) 
may become important. Thus condition (3.32) should be treated with caution. For instance, when 
(/3,w) is (2, 0.99), we found c\ and c 2 to be 1.6142 and 1.6355. In this case, uj 2 c 2 /c\ > 1. Thus 
according to (3.32), (fi,f 2 ) and (fi,— f 2 ) can not be matched. But our numerical results show 
otherwise [2]. 

Theorems 1 and 2 can also be used to study the matching of dark solitons which exist in coupled 
nonlinear Schrodinger equations (defocusing case) [12]. In this case, our results on the matching 
of some classes of dark solitons indicate that such matchings are impossible since conditions (2.68) 
can not be satisfied. We suspect that any dark solitons can not match each other to form widely- 
separated dark-soliton trains. 



4 Discussion 



The results in this paper can be readily applied to general nonlinear wave systems for the construc- 
tion of widely-separated multiple permanent-wave trains. Such wavetrains geometrically look like 
a superposition of individual permanent waves. This is somewhat analogous to the superposition 
principle of solutions in a linear system. But the difference here is that, due to the nonlinear na- 
ture of Eq. (2.2), those individual permanent waves have to be properly spaced (according to Eq. 
(2.21) or (2.68)) in order to form a wavetrain. When such wavetrains exist, one important question 
is their stability. For the coupled nonlinear Schrodinger equations, we indicated in [2] that they 
are all unstable. For certain Ginzburg-Landau and coupled-nonlinear-Schrodinger type systems, 
Malomed argued that multi-pulse trains exist and are stable by an approximate method based on 
the variational principle and effective potential ([13], [14]). Such results need to be viewed with 
caution due to the approximations involved. The clear evidence that some multi-pulse waves are 
stable can be found in the experimental results on binary fluid convection ([15]) and the numerical 
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results on subcritical Ginzburg-Landau equations ([16]). We will investigate those systems in the 
near future. 
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(b) 

Figure 1: Solitary waves matched by two wave and daughter 
wave solutions. The solid curves are ri(x), and the dashed 
curves are r2(x). In (a), (/?, w) = (2/3, 0.85) which is close to 
the curve (3.29) with m = 0; in (b), ((3,u>) = (2,0.6) which 
is close to the curve (3.29) with m = 1. 
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Figure 2: A solitary wave formed by mixed matching be- 
tween a wave and daughter wave solution and another soli- 
tary wave of different type. Here ((3,w) = (2/3,0.78), close 
to curve (3.29) with m = 0. The solid curve is r\(x), and 
the dashed curve is r2(x). 
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